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ABSTRACT 


This  report  develops  a  multidimensional,  dynamic  analysis  of  solid  state 
avalanche  diodes.  Well-established  electromagnetic  concepts  are  applied 
to  a  widely  used  model  of  the  diode  and  reveal  a  discrete  spectrum  of  new 
small-signal  modes.  The  approach  used  enlarges  the  conventional  per¬ 
spective  and  has  permitted  the  discovery  that  at  least  one  of  these  new 
modes  appears  to  possess  a  high  power  capability  (associated  with  its 
two-terminal  negative  resistance)  which  has  been  partially  realized  ex¬ 
perimentally.  The  lowest-order  mode  contains  all  the  results  of  prior 
quasi-static  theories  on  the  normal  IMPATT  mode,  plus  additional  infor¬ 
mation  which  is  used  to  delineate  the  range  of  validity  of  the  quasi -static 
results.  Formal  discrepancies  are  uncovered  between  the  usual  quasi¬ 
static,  one-dimensional  result  for  diode  impedance  used  in  solid  state 
studies  and  the  dynamic  multidimensional  result  for  the  normal  IMPATT 
mode  developed  from  microwave  circuit  theory.  However,  these  discrep¬ 
ancies  are  numerically  quite  small  except  in  certain  narrow  frequency 
bands. 
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I.  INTRODUCTION 


A  THEORY  OF  MULTIPLE  MODES 
IN  AVALANCHE  DIODES 


This  report  develops  a  multidimensional,  dynamic,  small-signal  analysis  of  solid  state 
1 

avalanche  diodes  by  applying  well-established  electromagnetic  concepts  to  a  widely  used  ava¬ 
lanche  diode  model.  Prior  theories  are  one-dimensional  and  quasi-static,  and  take  into  account 
only  longitudinal  (z)  variations  of  the  fields  and  currents  within  the  semiconductor.  The  result 
is  a  single  quasi-static  avalanche  mode.  In  the  present  work  radial  (r)  and  z -variations  are 
considered,  and  the  result  is  the  existence  of  a  spectrum  of  avalanche  modes  in  the  semicon¬ 
ductor.  The  lowest-order  mode  is  a  "volume"  mode  (i.e.,  a  mode  in  which  the  fields  and  cur¬ 
rents  do  not  vary  significantly  throughout  any  diode  cross  section)  and  is  just  that  mode  revealed, 
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in  part,  by  one-dimensional  theory.  The  remaining  modes  are  apparently  surface-wave  modes 
(i.e.,  modes  in  which  the  field  and  current  amplitudes  inside  the  semiconductor  decrease  signif¬ 
icantly  with  distance  from  the  diode-air  surface),  or  similar  to  surface  waves.  The  idealized 
solid  diode  configuration  with  its  contacts  may  be  pictured  somewhat  as  a  parallel-plate  capacitor 
or  a  radial  waveguide  with  perfectly  conducting  walls,  each  completely  filled  with  a  cylinder  of 
semiconductor  material.  The  dimensions  of  the  diode,  however,  are  several  orders  of  magni¬ 
tude  smaller  than  ordinary  capacitors  or  radial  waveguides  (see  Fig.  1). 


x  *-L/2 


*  *  L/2 


$  ■  0 


a  *  CONST 
0  «  CONST 


a  ■  /3  *  0 

•HOLE  CURRENT - 


■  ELECTRON  CURRENT  • 


x  *  b 

I 

I 


■  qPVp 
-qnV„ 


DRIFT 

ZONE 

HI 


AVALANCHE 

ZONE 

r 


DRIFT 

ZONE 

H 


Fig.  1.  Multiple  uniform  layer  model  of  avalanche  diode.  The  crosshatched 
sections  bordering  the  active  regions  represent  the  highly  conductive  bulk 
material. 

In  device  work,  mode  has  usually  signified  a  distinctive  pattern  of  device  behavior.  In  this 
report,  avalanche  mode  is  used  to  signify  a  self-consistent  field  and  current  distribution  within 
the  space-charge  region  of  the  diode  which  is  consistent  with  the  prevailing  boundary  conditions 
at  the  planar  interfaces  between  the  space-charge  region  and  the  remaining  material  of  the  diode. 
Each  such  mode  in  the  preceding  sense  turns  out  to  be  connected  with  a  distinctive  pattern  of 
device  behavior. 

The  avalanche  modes  are  obtained  by  deriving  and  solving  the  fundamental  small-signal  field 
equation  for  avalanching  carrier  interactions  with  electromagnetic  fields.  This  turns  out  to  be 
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Fig.  2.  Simplified  model  of  one  type 
of  packaged  diode. 


a  fourth-order  wave  equation  for  which  the  introduction  of  a  spatially  and  temporally  dispersive 
equivalent  conductivity  greatly  simplifies  the  method  of  solution.  The  avalanche  modes  within 
the  diode  are  coupled  to  the  purely  electromagnetic  modes  in  the  semiconductor  package  and 

mount  (see  Fig.  2).  For  many  package  and  mount  geometries,  these  modes  can  be  approximated 
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by  a  TEM  radial  waveguide  mode,  with  all  higher-order  radial  waveguide  modes  cut  off.  The 
modes  in  the  package  and  mount  are,  in  turn,  coupled  to  the  normal  electromagnetic  modes  of 
the  waveguide  or  transmission  line  containing  the  packaged  diode. 

The  approach  used  in  the  present  report5  enlarges  the  conventional  perspective  on  solid 
state  avalanche  diodes  in  a  manner  which  permits  rigorous  consideration  of  complex  circum¬ 
stances  which  are  not  directly  accessible  through  present  one-dimensional  theory.  From  this 
perspective  one  may  investigate  the  influence  on  avalanche  diode  operation  of: 

(a)  Higher-order  avalanche  modes  (which  apparently  possess  high  power 
generation  capability  at  microwave  frequencies), 

(b)  Noncircular  diode  cross  sections, 

(c)  Doughnut  or  ring  shapes^  (which  are  used  to  reduce  the  diode’s  thermal 
resistance), 

(d)  Nonuniform  radial  temperature  and  as  well  as  current  profiles  induced 
by  thermal  effects,*7 

as  well  as 

(e)  Determine  the  diode-package-mount  interaction  taking  into  account  the 
cut-off  higher-order  radial  waveguide  modes, 

(f)  Delineate  the  range  of  validity  of  the  conventional  one-dimensional,  quasi¬ 
static  theory  of  the  normal  IMPATT  mode,  and 

(g)  Analyze  potentially  self-resonant  IMPATT  diodes  (which  are  analogous  to 
self -resonant  LSA  diodes8). 

The  higher -order  avalanche  modes  appear  for  the  most  part  to  be  similar  to  surface  waves 
and  are  referred  to  as  SWIMPATT  (for  surface-wave  IMPATT)  modes,  and  are  discussed  at 

Q 

length  elsewhere.  The  modes  are  ordered,  as  is  the  case  of  waveguide  theory,  from  smallest 
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(lowest-order)  wavenumber  to  larger  (higher-order)  wavenumbers.  The  SWIMPATT  modes  are 

basically  a  two-dimensional  avalanche  phenomena.^  They  are  somewhat  analogous  to  the  recently 
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discovered  two-dimensional  wave  effects  in  Gunn-effect  diodes.  The  analogy  is  imperfect, 

however,  because  (1)  the  above-mentioned  Gunn-effect  diodes  are  presently  designed  in  the  co- 
planar  (or  sheet  diode)  configuration,  while  the  avalanche  diodes  under  discussion  are  assumed 
to  be  in  the  sandwich  configuration  described  at  the  beginning  of  this  section,  and  (2)  the  Gunn- 
effect  diodes  are  believed  to  operate  in  their  lowest-order  mode  while  SWIMPATT  are  higher- 
order  modes.  The  present  report  will  develop  the  fundamental  theory.  Topics  (b),  (c),  (e),  and 
(g)  will  be  discussed  in  Ref.  14  as  an  application  of  the  theory  developed  here.  Topics  (a)  and  (d) 
will  be  developed  in  Refs.  9  and  15.  The  present  report  will  consider  topic  (f)  and,  in  outline 
form,  topic  (a). 


II.  MODEL  AND  BASIC  EQUATIONS 


The  space-charge  region  of  the  avalanche  diode  is  modeled  as  multiple,  uniform,  tandem 


layers 


1,16 


(shown  in  Fig.  1)  within  which  a  ,  the  electron  ionization  rate  (which  is  assumed  for 


simplicity  equal  to  cv^,  the  hole  ionization  rate)  is  assumed  to  be  independent  of  the  longitudinal 
position  z.  It  is  also  assumed  that  the  material  properties,  the  DC  electric  field  Eq,  and  the 
DC  electron  and  hole  longitudinal  currents,  JnQ  and  Jpo,  are  independent  of  the  transverse  co¬ 
ordinates  within  the  diode.  The  space-charge  region  is  assumed  to  be  terminated  on  both  sides 
by  highly  conductive  material. 

In  the  avalanche  zone,  a  =  a  -  a  is  a  non-zero  quantity,  while  the  drift  zones  are  modeled 

n  P  _  A 

with  a  -  0  (see  Fig.  1).  The  electron  and  hole  velocities,  v  =  v  z  and  v  =  — v  z  ,  are  taken 
B  n  n  o  p  p  p 

to  be  their  scattering-limited  values  (where  v^  and  v^  are  positive  constants,  and  zq  is  a  unit 

vector  in  the  positive  z-direction).  Although  v  ^  v  ,  and  a  ^  a  in  Si  and  Ge  diodes,  it  is 

116  ^  P  ^  P 

assumed  for  simplicity  '  that  vn  =  v^  =  v,  and  o?n  =  =  a  (which  holds  for  GaAs  diodes). 

The  continuity  equations  for  the  electron  and  holes  are 


d  n  1  y  J  , 

9F  -  q  v  •  Jn  +  8 


(1) 


dp  _  -1 

at  q 


V 


J  +  g 
p  6 


(2) 


where  g  =  av  (n  +  p)  is  the  generation  rate,  =  —  qv^nz^  and  =  — qv^pz^  are  the  electron  and 
hole  current  densities,  q  is  the  electronic  charge,  and  n  and  p  are  the  electron  and  hole  den¬ 
sities.  The  effects  of  ohmic  radial  currents  are  considered  in  Ref.  9. 

The  small-signal  time-harmonic  approximation  for  the  electric  field  is 


e  =  E  +  E  e^wt  (3) 

o  > 

where  Eq  =  Eqzo  is  the  DC  component  of  e,  E  is  the  RF  component  of  e  with  a  frequency  cj 
radians  per  second,  and  |E  |  »  |E|.  A  similar  description  applies  to  the  electron  density 
N  =  nQ  +  n  elcJt,  and  the  hole  density  P  =  pQ  +  p  e^wt.  The  ionization  rate  is  approximated  by  the 
truncated  Taylor  series 

a(e)  &  aQ  +  «qEz  (4) 


t  This  differs  from  the  multidimensional  avalanche  phenomena  in  spherical  diodes  which  have 
recently  been  discussed.^ 
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expanded  about  Eq  with  E^  =  E  •  zq,  and  a Q,  =  dce/de,  evaluated  at  e  =  Eq.  Equations  (1) 

and  (2)  may  then  be  used  to  obtain  the  continuity  equations  for  the  RF  carrier  densities 


D  n  —  a  v  —  cE  =  0 
n  op  z 


D  p  -  a  v  -  cE  =  0 
p  on  z 

where  only  first-order  RF  terms  are  retained,  c  = 

n  8t  Bz  o 


~ao(Jno  +  Jpo)/q’  and 


(5) 

(6) 

(7) 


~  a  a 

d  =  ~  -  v  TT  -  01  v 
p  at  az  o 

In  the  one -dimensional  problem,  it  is  sufficient  to  add  Poisson's  equation 

aE_ 


E  = 


■w  * (p_n) 


(8) 


(9) 


to  obtain  a  complete  system  of  three  equations  in  three  unknowns  n,  p,  and  E^.  In  the  multi¬ 
dimensional  problem,  this  system  is  not  complete  because  of  the  existence  of  additional  field 
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components.  In  general,  this  requires  the  full  Maxwell  field  equations,  which  are 


V  x  E  = 


8H 

at 


(10) 


V  x  H  =  Jn  +  Jp  +  €  ^ 


(H) 


where  the  constitutive  relations  D  =  eE  and  B  =  pqH  have  already  been  used.  Equations  (5),  (6), 
(10),  and  (11)  constitute  a  complete  set  of  equations.  Equations  (5)  and  (6)  may  be  combined  to 
give 


[DD  -~(o'v)]n=c[D  +  a  vl  E 

1  n  p  o  J  lp  oJz 


[D  D  —  (a  v)  ]  p  =  c  [D  +  a  v]  E 
lnp  '  o  '  J  ^  ln  oJz 


(12) 

(13) 


Multiplying  Eqs.  (12)  and  (13)  by  —  qv,  and  summing 


7  3E 

[D  D  —  (a  V)  ]  [  J  +  J  ]  -  Za'  J  v 
1  n  p  o  1  1  n  pJ  o  o  at 


(14) 


is  obtained. 

Equations  (10)  and  (11)  can  be  combined  to  obtain 


atfc 


2  2 


(15) 


Operating  on  Eq.  (15)  by  the  operator  —  a^v  to  eliminate  Jn  and  J^,  it  follows  that 


2 

^ _ 

lD„Dp  -  V  x  V  x  E  ,  -2po«iJov  —f-  eo  -  «po  [DnD  -  «,ov)Zl  2-f 

at  o  t 


(16) 


Equation  (16)  is  the  fundamental  wave  equation  for  the  avalanche  diode  model  considered  in  this 
report.  Generalizations  to  the  important  case  of  ohmic  and  diffusive  radial  currents  are  con¬ 
sidered  in  Ref.  9.  Equation  (16)  is  a  fourth-order  partial  differential  equation  and,  in  general. 
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difficult  to  solve  directly.  For  the  circumstances  being  considered,  it  is  useful  to  take  advantage 
of  the  fact  that  the  and  operators  are  constant  coefficient  and  commute  with  the  operator 
(V  x  V  x).  Consider  electric  fields  which  are  describable  by 


E. 


(17) 


where 


then 


j(cL)t-K.z) 

E.  =  S.(r ,(p)  e 


(18) 


[D  D  -  (a  v)2]  E.  =  F.(o>,  K.)  E.  (19) 

1  n  p  '  o  1  l  l'  l  l  ' 

where  F.(co,  K.)  =  [ j(co  —  vk.)  —  a0vl  [j(w  +  vK.)  ~  aQv].  the  above  results,  along  with 

0/8t  =  jw,  Eq.  (16)  becomes 


V  x  V  x  E.  =  f . 

1  o  1 


(20) 


where  the  equivalent  permittivity  matrix  €.  is  defined  by 


£  .  = 

1 


I'  *  —  1 

L  J 


(21) 


in  which  €  is  the  standard  dielectric  permittivity  of  the  material  and  a.  is  the  equivalent  conduc¬ 
tivity  defined  by 


cr.  = 

l 


2 


a; 


—  Zi c ova  '  J 
J _ o  o 

+  2jo;  oqv  —  (K.v) 


2 


(22) 


with  J  =  J  +  J  .  a.  is  the  conductivity  in  the  sense  that  it  can  be  shown  from  Eq.(14)  that 
o  no  po  i  u 


J  .  +  J  .  =  cr.E  . 
ni  pi  i  zi 


(23) 


Equation  (20)  may  be  written  directly  in  terms  of  an  equivalent  conductivity  matrix  <r.  as 


V  x  V  x  E.  =  —  ju>|jL  [5.  +  jw£l]  *  E. 


(24) 


where  1  is  the  unity  matrix  and 


ai 


ai  0  0 

0  0  0 

0  0  0 


(25) 
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The  influence  of  highly  conductive  regions  sandwiching  the  space -charge  region  is  approx¬ 
imated  by  using  the  boundary  conditions  appropriate  to  infinitely  conductive  material,  i.e.,  the 
radial  component  of  E  is  zero  at  z  =  ±L/2.  The  fields  are  expressible  in  terms  of  TE  and  TM 

cylindrical  wave  modes  which  are  similar  in  some  respects  to  the  TE  and  TM  radial  waveguide 
17 

modes.  In  addition,  the  TE  radial  wave  modes  are  uncoupled  to  the  avalanche  currents  and, 
hence,  are  the  same  as  for  a  "cold"  dielectric  cylinder  bounded  by  metal  caps.  Thus,  only  TM 
radial  wave  modes  (Hz  =  0)  will  be  considered.  In  Ref.  14,  the  effect  of  variations  of  the  fields 
and  currents  with  azimuthal  angle  is  considered  in  connection  with  the  influence  of  noncircular 
diode  cross  section. 

The  TM  radial  wave  avalanche  modes  have  non-zero  field  components  Ez,  H^,  and  E^  (when 

there  are  no  azimuthal  variations)  and  are  also  similar  in  many  respects  to  the  well-known  small  - 

1 8 

signal  TM  modes  in  an  axially  directed  electron  beam  of  finite  radius.  When  a  TM  radial  wave 
avalanche  mode  has  "sufficiently  small"  radial  variations  (and  no  azimuthal  variations),  then  the 
radial  electric  field  component  E^  will  be  shown  in  a  later  section  to  have  a  negligible  magnitude 
from  the  point  of  view  used  in  most  computations.  The  resulting  avalanche  mode  can  be  envi¬ 
sioned  to  be  essentially  a  TEM  (in  the  r-direction)  wave  with  only  Ez  and  H  ^  as  significant  non¬ 
zero  field  components.  This  quasi-TEM  mode  for  which  the  z-variations  are  derivable  from 
one-dimensional  theory  (which  conventionally  explicitly  discusses  only  the  Ez  field  component) 
is  the  normal  IMPATT  mode.  The  field  component  for  the  normal  IMPATT  mode  is  readily 
derivable  in  the  terms  of  the  z -component  of  the  total  current  density  J,p  =  Jn  +  Jp  +  ju>Ez  which 
is  shown  in  Sec.  V  to  be  essential  independent  of  r  and  z  (as  is  already  widely  known).^  The 

TEM  radial  wave  avalanche  mode  differs  most  significantly  in  mathematical  form  from  the  TEM 
17 

radial  waveguide  mode  (which  has  Ez  and  as  the  only  non-zero  field  components)  by  virtue 
of  the  z-variations  in  Ez  which  the  space-charge  waves  in  the  avalanche  diode  produce.  The 
field  component  is  independent  of  z  in  both  the  TEM  avalanche  mode  and  the  TEM  radial  wave¬ 
guide  mode.  By  contrast,  it  may  be  noted  that  the  quasi-TEM  radial  wave  mode  (with  negligible 
value  for  Er)  in  an  axially  directed  electron  beam  of  "large"  finite  radius  does  exhibit  z-variations 
in  Ez  (Ref.  18). 

Solution  of  the  associated  boundary  value  problem,  posed  below,  yields  a  discrete  spectrum 
of  modes  within  the  diode,  each  of  which  has  particular  field  and  current  distributions  (which 
vary  with  frequency)  within  the  diode.  The  relative  amplitudes  of  the  modes  would  be  determined 
by  imposing  boundary  conditions  at  the  "rim"  of  the  semiconducting  pill,  i.e.,  at  r  =  R,  in  the 
same  manner  as  is  done  for  radial  waveguides  to  reflect  the  diode -circuit  interaction. 

III.  AVALANCHE  ZONE 

From  the  basic  equation,  Eq.  (24),  it  follows  without  further  assumptions  that  TM  waves  of 
j(o)t-k.z ) 

the  form  £>z-(r)  e  1  must  satisfy 


t  Using  Ampere's  rule  in  the  form 

p27r  p27r  pR 

Jo  Wio  Jo  JTrd*dr 

one  obtains  =  rJrp/2  inside  the  diode  and  R^Jrp/2r  outside  the  diode  (but  nearby),  where  R 
is  the  diode  radius. 


6 


(26) 


j_  _a_ 

r  8r 


('£) 


+  T^6  =0 
a  c  z 


where  the  subscript  i  has  been  dropped  from  6  for  simplicity.  Note  that  the  identity  V  x  V  x 
E  =  V ( V  •  E)  —  V2E  is  not  required  in  obtaining  Eq.  (26),  and  further  does  not  reduce  the  alge¬ 
braic  complexity  since  V  •  E  ^  0.  The  remaining  symbols  are  defined  by 


T  2  (K02-Kj2)(4-K.2) 

<Km  +  Ka2-Ki2> 


(27) 


2  2 
K  =  to  u  c 
o  ro 


=  (w/v)Z  +  2jaQ(w/v)  -K2 


(28) 

(29) 


Ka2,(o;a/v)2  .  u,a2  =  2«;JovA  .  (30) 
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K  is  the  standard  dielectric  wavenumber,  ±K  are  the  wavenumbers  calculated  by  Misawa 
o  m  J 

(along  with  zero)  for  one -dimensional  space-charge  waves  in  avalanching  carrier  streams,  and 

19 

a is  the  avalanche  frequency  originally  derived  by  Gilden  and  Hines  below  which  the  equiv¬ 
alent  diode  conductance  is  positive  in  their  theory. 

Equation  (27)  can  be  rewritten  as 


K4  -  K.2(K  2  + 
i  i  m 


Ko2  -  Ta2) 


+  [K2K2  -T2(K2  +  K2)] 
1  o  m  a  m  a  J 


(31) 


from  which  it  is  clear  that  there  are  four  values  of  K.  corresponding  to  T,  with 


K =  — K.  _ 
ll  i2 


K =  -K.  . 
i3  i4 


(32) 


For  simplicity  the  first  subscript,  i,  of  the  K’s  in  Eq.  (32)  will  be  dropped.  Thus,  the  only  non¬ 
singular  solution  that  includes  r  =  0  is^ 

i  “jK.z 

Ex=Jo(Tar)  X  Aie  1  <33) 

i=l 

where  J  (Tr)  is  the  zeroth  order  Bessel  function  of  the  first  kind,  A.  are  constants,  and  the  e^a  * 
o  1 

dependence  is  understood.  From  Eqs.(10)  and  (33),  it  follows  that 


and 


H  =  ja>£T  J,(T  r)  V 
J  a  1  a  '  Li 

i=l 


A.  e 


■jKjZ 


k2-k2 

O  1 


E  =  jT  J .  (T  r) 
r  J  a  1'  a 


K.A.  e 
i  i 


-jKjZ 


K  2  -  K  2 

O  1 


(34) 


(35) 


where  J^(Tr)  is  the  first-order  Bessel  function  of  the  first  kind.  Equations  (33)  to  (35)  describe 
the  only  non-zero  field  components. 

t  Ring  (or  doughnut)  diodes  require  both  Bessel  functions  of  the  first  and  second  kind  for  a  com¬ 
plete  description. 
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Equations  (12)  and  (13)  can  be  used  to  obtain 


-jK.z 

J„-W’  Z  1 

i=  1 


(36) 


J  =  J  (T  r)  7  a  .A.  e 
p  o  a  '  U  pi  i 


-jKjZ 


(37) 


i=  1 


where  o^.  =  (cr./2)  (1  +  vK./u),  and  c  .  =  (a./ 2)  (1  —  vK^/oj).  Note  <rn .  +  CTpj  =  CTj-  The  customary 
boundary  conditions  (these  are  not  altered  by  the  presence  of  drift  zones)  which  are 


J  =  0  at  z  -  —  — 
n  2 


J  =  0  at  z  =  +  — 
P  2 


are  symmetric  so  that  when  the  drift  zones  are  symmetric  about  z  =  0 
Jn(z)  =  Jp(-z) 

It  follows  from  this  symmetry  that 


Ai  =  A2 


A3=  A4 


Equations  (33)  to  (35)  then  simplify  to 

Ex  =  J0(Tar)  2A1  [cos(K1z)  +  M  cos(K3z)J 


(38) 

(39) 

(40) 

(41) 

(42) 


H  =  Zjwt  T  A.J,(T  r) 
c p  J  a  1  1'  a 


cos(K.z)  M  cos(K.z) 
+ 


K2  -K2 
o  1 


k2-k2 


(43) 


where 


Er  =  2jTaAlJl(Tar) 


sin(K^z)  MK^  sin(K^z) 


K  2  -  K  2 
o  1 


Ko  "K3 


M  =  7T  = 
Al 


u)  cos©^  +  jvK^  sinG^I  r  a;  -(vK^f  +  2jcjVQ'o 

u)  cos©-  +  jvK-  sin©.  I  2  /  xr  \2  t  -•  u 

3  J  3  3 J  l a;  — (vK^)  +  2ja;Va'o 


(44) 


(45) 


with  ©^  =  k^L/2,  and  ©^  =  K^L/2.  In  the  absence  of  symmetry,  the  results  are  more  complex. 
IV.  DRIFT  ZONE 

When  o'  and  a ’  are  almost  but  not  exactly  equal  to  zero,  Eqs.(5)  and  (6)  apply  in  the  same 
form  as  for  the  avalanche  zone.  Thus,  the  results  for  the  drift  zone  are 


-j/3.z 

Ex  =  Jo(Tdr)  Z  Gie  1 

i=  1 


(46) 
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(47) 


%  =  ^Vi(Tdrl 


E 


\(AYV 

t  /  piGi  \  -j^i2 

i=l'  o  pi  ' 


(48) 


where  =  ~@2’  ^3  =  ancl  are  constant  amplitudes.  However,  there  is  not  the  same  sym¬ 
metry  as  in  the  drift  region,  so  ^  and  G^  ^  G^.  The  equation  for  remains  in  the  same 
form,  i.e. 


Td  = 


<Ko-gi2)<Km-Pl2> 

<Km  +  Ka 


(49) 


But  now 


K  —  A  —0 
a  a 


because  a '  — •  0 
o 


K  1 


p2  +  A  .  -0 

m  m 


because  a  ,  a’  —  0 
o  o 


where  3  =  w/v  and  A  and  A  are  K  and  K  —  3  evaluated  at  arbitrarily  small  but  not  zero 
H  '  a  m  a  m  H  J 

values  of  a  and  a  ’  . 

o  o 

Thus,  Eq.  (49)  reduces  to 


(K2  -  p2)  (p2  -  IS2  +  A  ) 
T  2  _  o  i  i  m 

d  IS2  -  IS2  +  A  +  A 

i  m  a 


(50) 


Equation  (50)  can  be  rewritten  in  the  form 


t^]4  -  (^l2  (P2  +  K2  -T2)  +  pZ(K2  -T2)  =  0 


Ld '  H  '  o  d 

where  Aa  and  Am  have  been  taken  as  zero  for  simplicity.  Continuity  of  fields  at  x  =  ±L/2  re 
quires  that  T&  =  T^.  The  solutions  to  Eq.  (51)  are  then 


(51) 


P{  = 


(i  =  1,2,  3,4) 


(52) 


The  above  description  of  the  drift  zone  is  not  so  complex  as  that  required  for  the  avalanche 
zone.  For  the  unsymmetric  case,  there  are  four  unknown  coefficients  (the  G.’s)  for  each  of  two 
drift  zones,  plus  four  unknown  coefficients  (the  A.’s)  for  the  avalanche  zone  totaling  twelve  un¬ 
known  coefficients.  There  are  twelve  boundary  conditions  that  must  be  satisfied  at  the  various 


interfaces.  They  are  E  =  J 
J  r  p 

and  J  at  z 
P 


0  at  x  =  b,  E^  =  Jn  =  0  at  x  =  —a,  and  continuity  of  E^, 


H  .  J  , 
(p  n' 


±L/2.  This  information  is  necessary  and  sufficient  for  the  determination  of  the 
ratios  of  eleven  of  the  coefficients  to  the  twelfth,  plus  the  radial  wavenumber  T&.  Once  Ta  is 

determined,  then  the  longitudinal  wavenumbers  K.  and  (3.  can  be  calculated. 

2  2  2  1  1 

In  the  special  case  when  K  —  T&  =  (3  (which  corresponds  to  surface  waves  because  here 

|T  |  R  >  1  typically),  we  see  from  Eq.  (52)  that  a  degenerate  situation  occurs  in  which  case  eigen - 
a 

functions  A.  e  are  no  longer  valid.  To  analyze  this  special  case,  one  must  then  return  to 
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the  fundamental  avalanche  equation,  Eq.(16),  and  rederive  its  solutions  anew.  An  alternative 
approach  is  to  analyze  the  drift  zone  degenerate  case  with  =  0  to  begin  with,  as  opposed 

to  calculating  the  solutions  from  the  avalanche  zone  results  by  letting  q?q,  a'Q  —  0  as  we  have 
done  here.  However,  no  further  insight  is  gained  from  such  an  approach.  The  basic  wave  prop¬ 
erties  of  the  nondegenerate  surface  modes  are  outlined  in  Refs.  3  and  9. 

V.  THE  QUASI-TEM  RADIAL  WAVE  MODE  (NORMAL  IMPATT) 

In  this  section  it  is  demonstrated  that,  when  |  T&  |  r  «  1  for  0  <  r  $  R,  the  preceding  TM 

waves  reduce  to  a  quasi -TEM  cylindrical  bulk  wave  for  which  the  axial  component  of  electric 

field  is  that  derived  by  quasi-static,  one-dimensional  theory  for  the  normal  IMPATT  mode. 

The  first  part  of  this  demonstration  consists  of  deriving  approximate  expressions  for  the 

K.  's.  This  is  accomplished  by  solving  Eq.  (25),  for  K.  in  terms  of  T  ,  with  the  assumption  that 

|T  f  «  | K  |  .  For  typical  parameter  values  of  an  avalanche  diode  operating  at  microwave 

frequencies  (K  /K  )2  is  of  order  unity  so  that  |  T  | 2  1  +  (K  /K  )2  |  «  |  K  |  2  and  K  2  «  |  K  |  2. 

sl  m  sl  el  m  m  o  m 

The  solutions  to  Eq.  (31)  are  then 

Kl2  -  Kf  “Ko  -Ta  t1  +  <Ka/Km>2]  <53> 

K  2  =  K4  “  Km  -  Ko  +  Ta  t1  +  <Ka/Km)2)  <54> 

2  2  2 

Note  from  the  above  that  =  Km#  The  longitudinal  length  of  the  space-charge  region  is 

of  the  order  of  10  pm  or,  typically,  much  less.  Thus  cos(K^z)  ^  1,  sinfK^z)  ^K^z.  In  addition 
the  assumption  | T  |  r«  1  leads  to  J  (T  r)  =*  1,  J.(T  r)  =*T  r/2.  Then  Eqs.  (42)  through  (44) 

a  O  al  la  a 

reduce  to 


E 

z 

“  2A1 

[1  +  Mm  cos  (Kmz)] 

(55) 

jwer  , 

/  To 

2 

\ 

H 

<P 

A 

\  — 

J  [j^€rMm  cos(Kmz)] 

Ai 

1  *  <  W* 

\  K 
'  m 

-  A1(°'0  +  j“c)  r 

(56) 

/ 

'T  \ 

E 

r 

-  -jA 

,  sin  (K  z)  M  ( 
3  '  m  '  m  \ 

.O 

(Tar> 

(57) 

aQ  represents  a .  with  K.  =  0.  E z  is  now  in  the  TEM  form,  and  H^  is  in  the  TEM  form  except  for 
a  correction  term  which  is  negligible  for  |Ta/Km|  «  1.  E^  is  effectively  zero  because  it  is 
negligible  for  |  T  /K  |  «  1,  |  T  |  r  «  1. 

Equation  (45)  reduces  to 

A3  % 

Mm  =  AT  e(ja )  cose  —  VK  sine  ) 

1  J  mm  m 

where  0  =  K  L/ 2. 

m  m 

The  total  current  in  the  z -direction,  =  jcjcEz  +  can  be  obtained  from  Eqs.  (36), 

(37),  and  (42).  The  result  is 

JT  =  2A1JQ(Tar)  [ (cr1  +  jcoc )  cos  (K^z)  +  Mm(a3  +  Jcjc)  cos  (K3z)]  (59) 
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where  a  ^  =  a.  when  i  =  1,  and  a  ^  -  a.  when  i  =  3.  When  Eqs.  (53)  and  (54)  hold,  then  a  ^ 
and  (t^  Thus,  when  |  T&  |  r«  1,  then 

JT  =*  2A1(°r0  +  )  (60) 

Observe  that,  although  Jrp  is  independent  of  position,  it  is  a  function  of  frequency,  ionization 
rate,  saturated  carrier  velocity,  and  the  DC, current  density. 

If  one  attempts  to  approximate  the  preceding  quasi-TEM  radial  wave  by  a  cylindrical  wave 
which  is  TEM  without  approximation,  i.e.,  where  Ez  ^0,  and  ^  0  but  E^  =  0  (which  is  not 

possible  because  the  boundary  conditions  could  not  be  satisfied),  then  a  somewhat  unlikely  result 

20 

is  obtained.  It  can  be  shown  that  for  the  exact  TEM  radial  wave  the  unity  term  in  the  R.H.S. 
of  Eq.  (55)  is  associated  with  T&  given  by 

Ta  =  To  s  -^o(%  +  jwe)  Ko  <61> 

while  Ta  =  0  for  the  remaining  cosine  term,  i.e., 

E7  «  2A,  [J  (T  r)  +  M  cos  (K  z)]  (62) 

z  l  o  o  m  m 

Inasmuch  as  a  plays  the  role,  in  several  instances,  of  somewhat  of  an  effective  over-all  con¬ 
ductivity  (independent  of  K.)  for  the  quasi-TEM  (normal  IMPATT)  mode,  it  appears  plausible  to 
2  1 

take  Ta  as  given  by  Eq.  (61).  Then  Eq.  (55)  could  be  written  more  descriptively  as 

E  =  2A.J  (T  r)  [1  +  M  cos  (K  z)] 
z  loo'1  m  m 

Several  important  properties  of  Tq  are  readily  determined. 

(a)  If  cl  »  a;a  (and  typically  >  2a  qv), 

then  T  2  «  K  2  =  real, 
o  o 

(b)  If  a  -  0, 

(K  a;  \  /  K  \  lix  € 

23^7)  <1+  (2S7)  J^r 

(c)  If  w  =  +  (2oov)2, 

then  Re(T2)  =  0. 
o 

The  quasi-TEM  radial  wave  mode  within  the  IMPATT  diode,  as  explicitly  displayed  by 
Eqs.  (55),  (56),  and  (57)  can  be  now  readily  contrasted  with  the  TEM  radial  electromagnetic  wave¬ 
guide  mode  in  the  diode  package  and  mount  which  is,  in  general,  described  by 


E  =  2AJ  (K  r) 

X  o  o 

(63) 

H  =  -K  2AJ,  (K  r)/ju>(i 
< p  o  lv  o  '  J  'o 

(64) 

2  2 

where  K  =  cl  u  e  .  This  reduces  to 

0^0  0 
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(65) 


E  =  2A 
x 

=  jcoeAr  (66) 

for  the  usual  case  where  K  r  «  1. 

o 

The  influence  of  the  IMPATT  diode's  finite  radius  on  its  impedance  is  determined  below 
to  be 


Zc  TaRJo<TaR) 
Zo  =  2Jl(TaR> 


(67) 


where  Z  and  Z  are  the  impedances  derived  from  the  one  and  multidimensional  theories,  R  is 
o  c  2 

the  radius  of  the  diode,  and  one  may  use  the  simple  estimate  of  T&  provided  by  Eq.  (61).  The 

results,  shown  in  Figs.  3  through  10,  demonstrate  that  significant  deviations  (ranging  from  5  to 
50  percent  or  higher)  from  the  quasi-static,  one-dimensional  results  occur  only  within  a  100-MHz 
band  surrounding  the  threshold  frequencies  of  negative  impedance  and  the  reactive  resonance  fre¬ 
quencies.  The  substantial  deviations  illustrated  may  be  easily  obscured  in  practice  by  dissipative 
losses  in  the  contacts,  package,  and  mount,  as  well  as  the  usual  tolerances  in  measurements. 


Fig.  3.  Relative  percentage  difference  between 
normal  IMPATT  resistance  (in  a  100-MHz  band 
surrounding  the  negative  resistance  threshold) 
as  calculated  from  conventional  quasi -static 
one -dimensional  and  dynamic  multidimensional 
theories.  Parameters  are  JQ  =  3760  amp/cm^, 
L  =  0.50  pm,  b  — L/2  =  0.92  pm,  and  v  =  8.5  X 
10°  cm/sec. 


Fig.  4.  Relative  percentage  difference  between 
normal  IMPATT  reactance  (in  a  100-MHz  fre¬ 
quency  band  surrounding  reactive  resonance) 
as  calculated  from  one -dimensional  quasi¬ 
static  and  dynamic  multidimensional  theories. 
Parameters  same  as  for  Fig.  3. 
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I00(p«rc«nt)  [  ,00  (percent) 


Fig.  5.  Same  type  of  plot  as  in  Fig.  3, 
using  parameters  JQ  =  10, 350  amp/ cm2, 
L  =  0.22  pm,  b  —  L/2  =  0.39  pm,  and  v  = 
8.  5  X  10&  cm/sec. 


Fig.  6.  Same  type  of  plot  as  in  Fig.  4, 
using  same  parameters  as  for  Fig.  5. 


Fig.  7.  Same  type  of  plot  as  in  Fig.  3, 
using  parameters  JQ  =  14,000  amp/cm, 
L  =  0.17  pm,  b  —  L/2  =  0.30  pm,  and 
v  -  8.5  x  10 “6  cm/sec. 


Fig.  8.  Same  type  of  plot  as  in  Fig.  4, 
using  same  parameters  as  for  Fig.  7. 
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Fig.  9.  Same  type  of  plot  as  in  Fig.  3, 
using  parameters  JQ  =  23, 530  amp/cm^, 
L  -  0. 12  urn,  b  —  L/2  =  0.21  jam,  and  v  = 
8. 5  X  10®  cm/sec. 


FREQUENCY  (GHz) 


Fig.  10.  Same  type  of  plot  as  in  Fig.  4, 
using  same  parameters  as  for  Fig.  9. 


VI.  IMPEDANCE  OF  THE  AVALANCHE  REGION 


In  this  section,  Eq.  (67),  which  explicitly  displays  the  dependence  of  the  impedance  of  a  diode 
operating  in  the  normal  IMP  ATT  mode  on  diode  radius  R,  is  derived.  We  find  that  the  usual 
quasi-static,  one-dimensional  result  for  impedance  used  in  diode  work,  Eq.  (7 1)  below,  differs 
numerically  over  narrow  frequency  bands  from  the  dynamic,  multidimensional  result  for  the 
same  quantity. 

For  convenience,  we  repeat  here  the  appropriate  equations  for  the  normal  IMPATT  mode 
(where  cos^Z)  =*  1,  and  =*  Km)  explicitly  including  the  radial  dependence.  These  are,  from 
Eqs.  (42 ),  (43),  and  (44), 


Ez  “2AlJo(Tar)[l  +Mm  cos(Kmz)] 

%  “2AlJl<Tar)  <%  +  i"e>/Ta 

E  =*  0 
r 

The  impedance  of  a  diode  operating  in  a  normal  IMPATT  mode  is  given  in  general  as 


(68) 

(69) 

(70) 


(71) 


where  1^  is  the  total  current  through  the  diode  of  radius  R  and 
.L/2 

L 

-L/2  z 


rL 

V  =  \  E_dz 

J-I 


(72) 


is  the  voltage  across  the  diode.  The  voltage  may  be  written,  from  Eqs.  (68)  and  (72),  as 


V  =  Jo(Tar)  Vo 


(73) 
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where 


V 


•f 


L/2 

L/2 


2A .  [1  +  M  cos  (K  z)]  dz 
1  m  m 


(74) 


is  the  voltage  arising  in  the  conventional  one-dimensional  theory.  It  is  clear  from  Eq.  (7  3)  that 
the  voltage  remains  an  ambiguous  quantity  until  an  appropriate  value  of  r  is  chosen.  The  total 
current  is,  from  Eq.  (59), 


■.  ■  ft  j« 


rd  (p  dr 


(75) 


Eqs.(71),  (73),  and  (75)  yield  the  conventional  one-dimensional  impedance  Zq  =  V^/ttR  icto  +  ju>e), 
when  J  (Tr)  »  \t  and  J  is  given  by  Eq.  (60). 
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The  impedance  of  a  one-port  (two-terminal)  is  given  in  microwave  circuit  theory  by  * 


7  -  _ L_ 

c'  in2 


where 


E  X  H*  •  dA 


(76) 


(77) 


which  represents  the  complex  electromagnetic  power  flow  through  the  enveloping  surface  A c 
(which  in  the  present  case  is  the  "exposed"  cylindrical  surface  comprising  the  diode-air  inter¬ 
face  at  r  =  R).  I  in  Eq.  (76)  represents  the  terminal  current  which  may  be  chosen  as  the  1^  used 
in  the  present  report.  The  total  current  is  related  to  the  magnetic  field  by  Ampere's  Law  [which 
is  simply  the  usual  integral  form  of  Eq.(ll)] 


I  =  £  H 
t  J  (p 


R  d  <p 


(78) 


where  azimuthal  symmetry  has  been  assumed,  R  is  the  outer  radius  of  the  diode  and  in 
Eq.(7  8)  is  evaluated  at  r  =  R.  Since  is  independent  of  <p  (for  circular  cross-section  diode), 
Eq.  (78)  becomes 


I.  =  2?rRH 
t  <p 


=  47rA1RJ1(TaR)  (oq  +  jw«)/Ta 
For  the  normal  IMP  ATT  mode,  Eq.(77)  becomes 


Pt  =  ^<Ezn;)Rd^dz 


(79) 


(80) 


where  the  integrand  of  Eq.  (80)  is  evaluated  at  r  =  R.  This  result  may  be  rewritten  as 

-b 


=  I*  V  (evaluated  at  r  =  R) 


(81) 
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where  Eqs.  (72)  and  (7  8)  have  been  used  along  with  knowledge,  from  Eqs.  (68)  and  (69),  that 
is  independent  of  z  and  Ez  is  independent  of  <p.  Equation  (76)  thus  leads  to 


Jo(TaR)  Vo 

Zc  ;  47rAiR2J1(TaR2)  (crQ  +  ja>f)/T£ 


(82) 


or 


Z 

c 


( 


ZJ^^R)  ) 


Z 

O 


[Eq.  (67)] 


It  is  not  unexpected  that  a  multidimensional,  dynamic  description  of  a  diode  should  yield  some 
sort  of  corrections  to  the  one-dimensional,  quasi-static  results.  Such  is  the  case,  for  example, 
for  the  classical  capacitor.  However,  such  a  theory  also  reveals  the  potential  existence  of  other 
modes,  as  mentioned  earlier,  which  can  never  be  obtained  from  the  one -dimensional,  quasi¬ 
static  theory. 


VII.  THE  DRIFT  ZONE  IN  THE  QUASI-STATIC  APPROXIMATION 


In  this  section  the  quasi-static  one-dimensional  description  of  the  drift  zone  interactions  are 
retrieved  from  the  general  theory  of  Sec.  IV. 

When  |  T  |  r  «  1,  0  <  r  ^  R 


J  =  V  a  .G.  e  1 
n  Lj  m  l 
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similar  to  the  avalanche  zone  where  B.  replaces  K.  in  a  .  and  a  ..  For 
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then  lj  »  (v/3.)  for  all  c j,  so  —  an2  -  cr  ^  -  <7p3  -  a  o  goes  to  zero  where  aQ  and  go  to 
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for  sufficiently  small  a  and  a’ .  From  this  it  can  be  shown  that  a  .  and  0  .  go  to  zero  while 
^  o  o  n4  p3  b 

crp4  and  an3  go  to  -jo>c. 
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The  boundary  condition  Jp  =  0,  at  z  =  b,  implies  from  the  above  that 
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range  of  values  for  z,  e  ^  1  with  considerable  accuracy.  Thus, 
reduce,  for  the  quasi-TEM  case,  to 
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=  0.  For  the  customary 
Eqs.  (46),  (47),  and  (48) 
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where  Gq  =  +  G^.  The  boundary  condition  Ef  =  0  at  z  =  b  implies  that  G^  =  G^  for  T  //3  suf¬ 

ficiently  small,  so  the  above  reduces  further  to  E  ^  G  +  G0  e~^z,  H  «  icjf  rG  /2,  and  E  0. 

J  x  o  3  <P  o'  r 

Thus  the  quasi-static  results  can  be  retrieved.  The  large  radius  result  of  Eq.  (67)  applies  equally 
well  to  the  drift  zone. 


VIII.  CONCLUSIONS 

This  report  has  presented  a  multidimensional  analysis  of  avalanche  diodes  which  exposes 

the  existence  of  new  modes  (TM  cylindrical  waves)  in  the  semiconductor  in  addition  to  that  ob- 
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tained  through  prior  one-dimensional  theories.  '  The  main  parameters  characterizing  the 

radial  variation  of  fields  and  currents  are  (1)  the  frequency,  (2)  the  avalanche  frequency,  as  de- 
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fined  by  Gilden  and  Hines,  and  (3)  the  wavenumbers  for  one-dimensional  space-charge  waves 
on  avalanching  carriers  as  determined  by  Misawa.*^  Associated  with  each  mode  are  four  plane 
waves  traveling  along  the  longitudinal  axis  of  the  semiconductor.  The  results  obtained  from  one¬ 
dimensional  quasi-static  analysis  turn  out  to  be  contained  in  a  quasi-TEM  cylindrical  "volume" 
wave  which  is  retrieved  from  the  lowest-order  mode  of  the  multidimensional  avalanche  mode 
spectrum.  The  total  current  for  the  quasi-TEM  mode,  which  in  prior  analyses  is  an  unknown 
constant,  has  negligible  spatial  variation  but  is  shown  to  be  a  function  of  the  DC  current  density, 
frequency,  ionization  coefficients,  and  the  saturated  drift  velocity.  The  power  transport  and 
the  impedance  have  been  calculated  from  E  and  H.  The  quasi-TEM  mode  is  the  normal  IMPATT 
mode. 

A  small  discrepancy  has  been  uncovered  between  the  quasi-static  results  for  the  normal 
IMPATT  mode  diode  impedance  conventionally  used  in  diode  research  [Eq.  (71)]  and  the  dynamic 
result  [ Eq .  (76)]  developed  from  multidimensional  theory.  Fortunately,  the  magnitude  of  the 
numerical  discrepancy  is  negligible  for  most  frequencies  of  interest  except  in  narrow  bands  of 
about  100 -MHz  width  where  it  may  be  obscured  in  practice  by  other  factors.  However,  the  ex¬ 
istence  of  the  discrepancy  does  leave  the  question  of  whether  its  magnitude  may  not  be  significant 
in  special  diode  applications  (e.g.,  at  submillimeter  frequencies)  or  for  diodes  of  other  types. 

The  mode  spectrum  has  been  determined  by  the  continuity  conditions  at  the  planar  bounda¬ 
ries  z  =  ±L/2,  and  the  boundary  conditions  at  the  planar  boundaries  z  =  +b  and  z  =  —a.  This 
corresponds  exactly  to  the  procedure  used  in  the  prior  one-dimerisional  theories.  As  in  the 
theory  of  radial  waveguides,  specific  conditions  at  the  "rim,"  i.e.,  at  r  =  R,  can  be  synthesized 
by  an  appropriate  superposition  of  avalanche  modes  assuming  mode  completeness  (which  has  not 
been  proved  here).  Generally  there  are  no  constraints  that  may  be  imposed  at  r  =  R  because 

this  will  vary  with  the  tuning  of  the  associated  microwave  circuit  and  particulars  of  the  local 
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structures.  An  exception  is  the  self-resonant  diode. 

The  spectrum  of  higher-order  modes  is  discrete  and  typically  involves  RF  field  and  RF  cur¬ 
rent  profiles  which  vary  strongly  with  radial  position.  It  is  reasonable  to  expect  some  coupling 
between  the  various  avalanche  modes  to  occur  as  the  signal  level  is  increased.  Thus,  a  rigorous 
large-signal  theory  of  avalanche  diodes  must  evaluate  the  effects  of  such  mode  coupling. 

Geometric  effects  and  higher-order  modes  are  considered  in  detail  in  the  next  two  parts  of 
this  work  in  Refs.  9  and  14.  However,  it  is  useful  to  describe  here,  in  a  qualitative  fashion,  the 
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unusual  but  practical  nature  of  a  higher-order  mode.  Consider  an  avalanche  diode  whose  cross- 
sectional  area  is  one  to  two  orders  of  magnitude  greater  than  those  of  conventional  avalanche 
diodes  operating  at  the  same  frequency.  The  capacitance  of  such  a  diode,  at  microwave  frequen¬ 
cies,  would  be  one  to  two  orders  of  magnitude  greater  than  the  permissible  level  for  normal 
IMPATT  operation.  However,  the  preceding  observation  assumed  a  uniform  RF  field  profile 
across  the  cross  section  of  the  diode.  If  instead  the  RF  field  has  an  exponential  taper  with  radial 
position,  as  is  approximately  true  in  a  SWIMPATT  mode,  then  the  effective  capacitance  would 
be  moderate.  Figure  11  illustrates  this  situation.  Further,  if  one  calculates  the  "effective  ac¬ 
tive  volume  for  the  SWIMPATT  mode,  taking  into  account  the  sharp  profile,  it  turns  out  to  be 
larger  for  this  surface  wave  than  for  the  corresponding  volume  wave  in  normal  IMPATT  diodes. 
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Fig.  11.  Radial  profiles  of  axially  directed  electric  field  for  (a)  normal-sized  normal 
IMPATT  mode  diodes,  and  (b)  large-area  SWIMPATT  mode  diodes. 


The  buildup  of  a  SWIMPATT  wave  may  be  conveniently  pictured  in  two  simplified  stages. 

The  first  stage  may  occur  in  a  time  period  of,  say,  10  nanoseconds  (for  an  X-band  diode)  during 
which  the  effects  of  diffusion  are  too  slow  to  be  noticeable.  In  this  stage,  the  fields  and  current 
quickly  build  up  to  a  state  in  which  they  have  a  steep,  nearly  exponential  variation  with  radial 
position,  from  their  peak  values  at  the  outer  diode  perimeter  to  almost  negligible  values  several 
tens  of  microns  further  inward  toward  the  center  of  the  diode.  The  power  output  from  such  a 
wave  is  of  the  same  order  as  that  from  the  corresponding  normal  IMPATT  diode  because  the 
effective  active  volume  of  the  semiconductor  is  of  the  same  order  in  both  cases. 

As  time  passes,  diffusion  occurs  because  of  the  steep  RF  radial  gradients  in  current  density. 
In  other  circumstances,  such  diffusion  would  be  small  since  the  carriers  would  be  swept  out  of 
the  thin  active  region  in  the  axial  direction  at  their  saturated  drift  speeds  (vg).  Their  radial 
excursions  would  be  very  limited  since  their  average  radial  velocities  due  to  diffusion  would  be 
very  much  smaller  than  v  .  However,  avalanching  carriers  continually  produce  electron -hole 
pairs  throughout  the  avalanche  zone.  Thus,  as  they  travel  along  their  relatively  small  radial 
excursions,  the  avalanching  carriers  would  produce  other  carriers  which  would  extend  the  radial 
excursion  and  in  turn  produce  other  carriers  to  further  extend  the  radial  inflow.  The  net  result 
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is  that  the  active  volume  increases  markedly  so  that  the  output  power  goes  up  at  least  an  order 

of  magnitude.  However,  the  inflow  of  carriers  increases  the  axial  electric  field  further  in  from 

the  perimeter,  which  raises  the  effective  capacitance  and  reduces  the  frequency  of  oscillation. 

The  picture  obtained  from  the  above  considerations  is  that  of  a  very  large  area  avalanche 

diode  which  produces  high  power  as  its  frequency  of  oscillation  sweeps  over  a  broad  range. 

Such  diodes  have  been  fabricated  in  the  past  and  their  reported  operation  was  considered  unex- 
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plained  in  numerous  details.  *  Now  some  details  may  be  accurately  explained  by  the  above 

description  and  by  further,  more  quantitative  information  to  be  developed  in  Ref.  9.  For  example, 

calculations  based  on  some  normal  IMPATT  thermal  considerations  yielded  estimates  of  their 

sweep  frequencies  from  two  to  three  orders  of  magnitude  off  from  the  observed  sweeps,  while 

22 

other  more  speculative  suggestions  yield  no  quantitative  estimates.  On  the  other  hand,  some 
simple  calculations,  presented  in  Ref.  9,  obtain  both  the  correct  order  of  magnitude  of  the  sweep 
frequencies  and  the  correct  shape  of  the  curve  of  sweep  frequencies  vs  time.  As  a  further  ex¬ 
ample,  the  existence  of  the  SW1MPATT  mode  and  its  characteristics  explain  the  descrepancy 
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between  Melick’s  work,  (which  demonstrated  that  the  RF  output  IMPATT  power  peaked  for  areas 

around  the  ’’normal  sized”  IMPATT  diode  and  decreased  significantly  with  increasing  diode  area), 

and  the  exceptionally  high  peak  power  obtained  from  the  extra  large  area  diodes.  These  large 

area  diodes  produced  one  to  two  orders  of  magnitude  more  RF  power  than  the  normal  sized 

IMPATT  diodes  (operated  at  about  the  same  frequency)  fabricated  by  the  same  workers  at  their 
24 

time  of  discovery.  The  multiplicity  of  pulsed  modes  experimentally  detected  at  the  same  time 
and  in  the  same  diode  by  Manasse  and  Shapiro25  is  also  consistent  with  the  multiple  mode  theory 
developed  in  this  report  and  Ref.  9. 

A  final  example  is  the  efficiency  predictions  of  DeLoach  based  on  skin  effects  in  the  diode 
substrate  which  indicate  a  marked  decrease  from  the  normal  efficiency  of  an  IMPATT  diode  as 
its  radius  exceeds  the  skin  depth  in  the  diode  substrate.^  Gilden  and  Moroney24  obtained  approx¬ 
imately  the  same  efficiencies  (5  percent)  for  36-mil-diameter  diodes  as  those  obtained  for  nor¬ 
mal  (4-mil)  IMPATT  diodes  by  the  same  workers  at  that  time.  This  apparent  inconsistency  is 
explainable  by  the  surface-wave  character  of  the  SWIMPATT  mode. 
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